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The discovery of topological insulators-*^ has 
sparked considerable interest in the study of topo- 
logical phases of matter. The hallmark of these 
novel phases is the emergence of topologically 
protected boundary phenomenzt^^, e.g. quantum 
pumping^ ^, surface states related to exotic mod- 
els from particle physics'* and quasi-particles with 
non-Abelian statistic^. Yet, realizations of these 
phases of matter are scarc^^^^. Additionally, 
one-dimensional systems without any additional 
symmetries are considered to lack such topologi- 
cal phenomena-*^^. In this Letter, we show, both 
theoretically and experimentally, that quasicrys- 
tals possess high-dimensional topological phases. 
We demonstrate that one-dimensional quasicrys- 
tals exhibit topological properties that were, thus 
far, thought to be limited to two-dimensional sys- 
tems. Using photonic quasicrystals, we observe 
localized boundary states, which manifest these 
topological properties. The topological nature of 
these boundary states is used to realize an adi- 
abatic optical pumping effect in which photons 
are uni-directionally pumped across the sample. 
Generalization to various types of quasicrystals 
and higher dimensions are also discussed, suggest- 
ing the existence of topological effects on surfaces 
of three-dimensional quasicrystals. This work 
opens a path to the realization of new types of 
topological phases, and provides new insights to 
the physics of quasicrystals. 

Topological phases of matter consist of various band 
insulators or superconductors that have gaps in their 
spectrunP. Two systems belong to the same topological 
phase if they can be continuously deformed from one into 
the other without closing energy gaps. Consequently, at 
the interface between two topologically distinct systems 
the energy gaps close by the appearance of boundary 
states which are localized on the interface. Classifica- 
tion of all the possible topological phases according to 
their symmetries and dimensionality was recently intro- 
ducecpS For example, in the absence of any symmetries, 
all the one-dimensional (ID) systems belong to the topo- 
logically trivial phase, while in two-dimensions (2D) there 
is a well known series of topological phases, namely those 
of the Integer Quantum Hall Effect (IQHEf^. 

Quasicrystals (QCs) are non-periodic structures that 
possess long-range ordei^i^. This quasi-periodic order can 
be seen as originating from periodic structures of a di- 
mension higher than the physical one. For example, the 
ID Fibonacci QC can be described as a projection of a 
2D lattice on a lin^^^. Remarkably, observed phenomena 



such as unconventional Bragg diffraction and the exis- 
tence of PhasonJiSHni Q^n be attributed to this higher 
dimension. Indeed, remnants of the higher dimension- 
ality appear as additional degrees of freedom (d.o.f.) in 
the form of shifts of the origin points of the quasi-periodic 
order. These d.o.f. discern between QCs with the same 
quasi-periodic order as they result in different patterns. 
However, they have no apparent influence on the bulk 
properties, and were therefore usually ignored. 

In this work, we show that due to these additional 
d.o.f. , QCs exhibit non-trivial topological properties 
that are usually attributed to systems of a higher di- 
mension. The topological properties of the QC manifest 
in two ways: (i) The existence of quantum phase transi- 
tions when continuously deforming between two topologi- 
cally distinct QCs; (ii) The appearance of robust bound- 
ary states which traverse the bulk gaps as a function 
of the aforementioned shifts. We demonstrate experi- 
mentally that ID photonic QCs exhibit such boundary 
states. Harnessing the topological nature of these bound- 
ary states, we generate a quantum pump in which pho- 
tons are transferred across the QC by shifting the pat- 
tern's origin adiabatically. 

Let us begin with a specific example of a QC, the ID 
Aubry-Andre (AA) modeP^ (also known as the Harper 
model). This is a ID tight-binding model in which the 
on-site potential is modulated in space. It is described 
by the Hamiltonian 

H{(p)tpn = t{lpn+l + V'n-l) + A COS {2TTbn + (j)) Ipn ■ (1) 

Here, is the wavefunction at site n, t is the hopping 
amplitude, A is the modulation amplitude of the on-site 
potential, and b controls the periodicity of the modula- 
tion. Whenever b is irrational the modulation is incom- 
mensurate with the lattice and the on-site term is quasi- 
periodic. Note that this model embeds the d.o.f. men- 
tioned above in the form of the phase (p of the spatial 
modulation. 

Fig. [T^ depicts the numerically calculated spectrum of 
the AA model as a function oi 4> ior t — 1, X — 0.5, 6 = 
(Vb + l) /2 (the Golden mean) and n = -49.. .49. Due 
to the incommensurate potential, the spectrum is broken 
into a fractal set of bands and gapJ^^. Note that as a 
function of the bands are almost unchanged, yet the 
gaps are crossed by a few modes. The states that reside 
within the gaps are boundary states, localized either on 
the left or on the right boundary of the system, as seen 
in insets (1) and (2), respectively, while the states within 
the bands are typically extended, as depicted in inset (3) . 
As we later show, these boundary states are the physical 
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FIG. 1: Observation of topological boundary states in an Aubry-Andre photonic quasicrystal. (a) The numerically calculated 
spectrum (marked by blue and red dots) of equation (jlj) as a function of the phase for t = I, X = 0.5, b = (\/5 + l) /2, and 
n — —49. ..49. The bulk of the spectrum remains fixed, whereas few modes, localized at the boundaries, sweep across the gaps. 
The insets depict the spatial density of typical eigenstates as a function of position along the ID lattice: (1) a left boundary 
state; (2) a right boundary state; and (3) an extended state within the bands, (b) An illustration of the experimental system 
- the photonic waveguide lattice. Light that is injected into one of the waveguides tunnels to neighboring waveguides as it 
propagates. The propagation is described by equation ([l]) where the propagation axis, z, takes over the role of time, (c)-(e) 
Experimental observation of the left boundary state for (p = n/2. Light was initially injected to a single waveguide (marked by 
red arrows). The measured outgoing intensity is plotted versus the position along the lattice, (c) An excitation at the middle of 
the lattice (site 0) results in a significant spread, (d) When the light was injected to the rightmost site (site 49), the excitation 
again shows comparable spread, (e) For an excitation at the leftmost site (site —49), the light remains tightly localized at the 
boundary, marking the existence of a boundary state. 



manifestation of the fact that the AA model belongs to 
a non-trivial topological phase. 

We implemented the AA model in an optical setup us- 
ing a quasi-periodic array of coupled waveguides^, as 
illustrated in Fig. [T|d (see Methods). In such systems 
light is injected into one of the waveguides, and is coher- 
ently coupled to neighboring waveguides as it propagates. 
The propagation of the light is described by an equation 
which is identical to a tight-binding model, where the 
propagation axis takes over the role of tim^^U By modu- 
lating the refraction index of the waveguides it is possible 
to directly realize the AA Hamiltonian of equation ([T]) . 

We fabricated an AA lattice with the above mentioned 
parameters and (j) = 7r/2. The experimental observations 
are depicted in Fig. [lj;-e. When the light was injected 
into a lattice site in the middle of the array, the outgo- 



ing intensity distribution showed a significant expansion, 
due to the overlap of the injected wavefunction with the 
extended bulk eigenstates. When light was injected into 
the rightmost lattice site it again showed considerable 
expansion, meaning that here too, the initial excitation 
overlapped mostly with extended bulk states. However, 
when the light was injected to the leftmost lattice site, 
the intensity distribution remained tightly localized at 
the boundary, with the maximum intensity found at the 
leftmost waveguide itself. This is a clear signature of the 
existence of a topological localized boundary state. 

A consequence of the topological nature of this model is 
that all boundary states which reside within the same gap 
belong to the same mode. This can be seen by following 
the eigenenergy of some boundary state as a function of 
(j). For example, take the right boundary state denoted 
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in inset (2) of Fig. [Tp,. Its mode is marked by red circles. 
This state remains localized on the right boundary as 
long as its energy remains within the gap. When the 
energy reaches the band, the state becomes extended. 
Notably, once the mode returns into the gap it appears 
localized at the left boundary. 

This property was used to realize adiabatic pumping of 
photons from one side of the lattice to the other. A con- 
venient platform for this feat is the "off-diagonal" version 
of the AA model which is described by the Hamiltonian 

Ho«{cj))i^n = t [1 + A COS (27r6n + cj))] i^n+i (2) 
+ < [1 + A cos (27r6(n - I) + (j))] tpn-i ■ 

While this model embeds its quasi-periodicity in the hop- 
ping term, it has similar topological characteristics as its 
previously discussed "diagonal" version [cf. equation Q]. 
The pumping takes place when is adiabatically swept 
along the propagation axis z. This is achieved by slowly 
modifying the spacing between the waveguides along the 
propagation axis, as illustrated in Fig. [2^. In our imple- 
mentation we chose t = 30/78, A = 0.6, b = (\/5 + l) /2, 
and n = 1...21, and the length of the sample to be 75mm, 
which are 19.5 tunneling lengths. Fig. [2]d depicts the 
spectrum of the system for these parameters as a function 
of (j). The pumping from the right boundary to the left 
one was carried out by sweeping (j) from 0.47r to I.Gtt. In 
order to observe different stages of the pumping process, 
we fabricated a set of 50 samples for which the light was 
allowed to propagate shorter distances within the modu- 
lation. Correspondingly, in the i"^ sample (p is modulated 
from OAtt to (0.4 -|- 1.2 • j/50)7r. For each sample, light 
was injected to the rightmost site and the output inten- 
sity distribution was measured. The collected results are 
summarized in Fig. [2j:. The obtained intensity distribu- 
tions are stacked incrementally according to their prop- 
agation distance, i.e. their final (j). Thus, we reconstruct 
the light's trajectory along the full adiabatic process. It 
is evident that the injected light was pumped adiabati- 
cally across the QC from one boundary to the othei'221. 

We now turn to establish our results theoretically. We 
start by showing that the observed boundary states are 
of topological origin by mapping the AA model to the 
lattice version of the 2D IQHE~. In the latter, elec- 
trons hop on a 2D rectangular lattice hopping ampli- 
tudes t and t' in the presence of a perpendicular mag- 
netic field, with b flux quanta threading each rectan- 
gle. Assuming one coordinate to be periodic and us- 
ing the Landau gauge for the magnetic field, the sys- 
tem can be described by the Hamiltonian 'HV'ri.fc = 
Xlfc [t{->pn+i,k + ipn-i,k) + 2t' cos{2'Kbn -f ka)Tl)n,k]y where 
k is the momentum along the periodic coordinate with 
lattice spacing a and n is the location in real space along 
the second coordinate. The energy spectrum of H is 
gapped, and each gap is associated with a quantized Hall 
conductance ajj = I'e^/h, with v an integer^ known as 
the Chern numbei'^. The inclusion of disorder and dis- 
tortions in the Hamiltonian does not alter an as long as 
the corresponding gap is maintained operP^HUl. Due to 



the fact that the energy gap must be closed in order for 
an to change its value, it can be used to classify different 
phases of the IQHE. Phases with different are said to 
be topologically distinct. 

The physical manifestation of a non-trivial topologi- 
cal phase (i.e. (Th ^ 0) is the emergence of robust chiral 
states along the edges of the sample. This general phe- 
nomena is shared by many topological phases, not only 
the IQHEPHH. in the IQHE, on each edge exactly |z/| 
edge states appear, and the energy of each edge state tra- 
verses the gap as k varies from —it /a to tt/o. The signs 
of the group velocity of these edge states are opposite on 
opposite edges. 

Turning back to the ID AA model [cf. equation ([T])], 
we shall now observe that it inherits its robust bound- 
ary states from the 2D IQHE. For A = 2t' the spectrum 
of H{(j)) can be viewed as the k^^ component of H at 
k = (j)/a. Therefore, by scanning from — tt to tt, the 
spectrum of H{cj)) reconstructs that of H. Consequently, 
the chiral edge states that traverse the gaps as a function 
of k appear now as ID boundary states that traverse the 
gaps with (j). Since these boundary states are of topolog- 
ical origin, the only way to eliminate them is to close the 
energy gap which they traverse. In particular, disorder 
which does not close an energy gap does not eliminate 
the corresponding boundary states. 

Note that the topology guarantees the existence of 
boundary states only for intervals of (j). Hence, it does 
not guarantee that for any QC pattern they indeed ap- 
pear. This can be seen, for example, in Fig. [2]3 where 
ioT (j) — one finds states localized on both boundaries, 
and for (f> = tt there are none. Additionally, changing the 
number of lattice sites or translating them (e.g. taking 
n = 5. ..25 instead of n = 1...21), may alter dramatically 
the intervals of for which boundary states appear. 

So far, in order to witness the topological nature of the 
AA model we had to scan cj). Ostensibly, one can asso- 
ciate Chern numbers only to the union of all the H((j)) 
Hamiltonians. Indeed, per H{(j)) a Chern density is as- 
signed, while the Chern number involves integration of 
the Chern density over all </>. However, we show that for 
QCs the Chern density is independent of <j)- Thus, the 
Chern number can be evaluated from the Chern density 
of any H{(j3). Since the same quantized Chern number 
is associated with any H{(j)), it topologically classifies it. 
This is somewhat analogous to the role of the Aharonov- 
Bohm flux in the IQHE. Formally it is integrated over, 
while in actual experiments there is no s uch fl ux and yet 
a quantized Hall conductance is obtained ^^ l ^^ l The asso- 
ciation of a Chern number with each QC is a key result 
of this work. 

Proof of the above statement appears in the Supple- 
mentary information-251. Here we show that the bulk spec- 
trum is also independent of (p. This simpler proof con- 
tains the essential ingredients of the former. Since H{(j)) 
has a band structure, the spectrum is insensitive to lattice 
translations in the thermodynamic limit. From equation 
(fTl) it is evident that translating the lattice by m sites 
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FIG. 2: Experimental observation of adiabatic pumping via topologically protected boundary states in a photonic quasicrystal. 
Here, we use the "off-diagonal" version of the Aubry-Andre model [see equation ([2|] where the quasi-periodicity is in the 
hopping term, (a) An illustration of the adiabatically modulated photonic quasicrystal, constructed by slowly varying the 
spacing between the waveguides along the propagation ajcis z. Consequently, the injected light experiences an adiabatically 
modulated Hamiltonian, _ffoff(0(z)), as it propagates and is pumped across the sample, (b) The spectrum (marked by blue and 
red dots) of the model as a function of the phase for t = 30/78, A — 0.6, 6 — (\/5 + l) /2, and n = 1...21. In the experiment (f> 
was scanned between OAn and I.Gtt as marked by the red dots. The insets depict the spatial density of a boundary eigenstate as 
function of the position at three different stages of the evolution: at <^ = OAtt the eigenstate is localized on the right boundary 
of the lattice. At = tt it is delocalized across the system, while at (p = I.Btt the state is again localized, but on the left 
boundary of the sample, (c) Experimental results: Light was injected to the rightmost waveguide (site 21) a,t z — {(j> ~ OAtt). 
The measured intensity distributions as a function of the position are presented at different stages of the adiabatic evolution, 
i.e. different propagation distances. It is evident that along the adiabatic evolution the light crossed the lattice from right to 
left. 



is equivalent to shifting (f) hy 2tt ■ {bm modi). Now the 
irrationality of b comes into play. For a rational b ~ p/q, 
{bm modi) has only q different values for all possible 
translations. Thus, the band structure is guaranteed to 
be invariant only for these q corresponding shifts of 0. 
On the other hand, for irrational 6, {bm modi) samples 
the entire [0, 1] interval and the bands are invariant for 
any shift of 0. 

Thus we arrive at the following conclusion: while in 
order to witness boundary effects the scanning over (j) is 
required, the topological indices can be associated with 
any instance of a quasi-periodic pattern, i.e. any given 
(j). These indices are, of course, the same for a given 
quasi-periodicity for all 0's. Thus the AA model is topo- 
logically classified. Consequently, two QCs with two dif- 
ferent 6's cannot be smoothly deformed from one to the 
other without closing the bulk gaps, since in the IQHE 
different b's result in different Chern numbers^^. 

Until now, wc focused on a specific model which we 
were able to map to the IQHE. However, our results could 
be easily generalized to any QC, such as the off-diagonal 
AA model and the Fibonacci QC. Moreover, our argu- 
ments apply to any dimension and any topological index 
without the need for establishing such a mapping^^. Con- 
sider a D-dimensional QC with a tight-binding Hamilto- 
nian with d quasi-periodic terms, either hopping or on- 
site. These terms result in d d.o.f. similar to the above (p. 
In the context of topological properties, these d.o.f. could 
be treated as extra dimensions, yielding an overall ef- 



fective dimension of D + d. Therefore the Hamiltonian 
may belong to a non-trivial D + d-dimensional topologi- 
cal class. In such a case, the boundaries of the physical 
QC host boundary states for intervals of the d.o.f. and 
a quantum phase transition takes place when deforming 
between QCs with different topological indices. 

In conclusion, we show that quasicrystals exhibit new 
types of topological phases, that were previously at- 
tributed only to systems of higher dimension. We demon- 
strate experimentally such topological effects in ID sys- 
tems using photonic quasicrystals. Moreover, we harness 
the topological nature of these ID quasicrystals to re- 
alize an adiabatic pumping effect, observed directly in 
the optical setup. The study of the novel topological 
phases in 2D and 3D may lead to the discovery of new 
surface phenomena in atomic and photonic quasicrystals, 
e.g. 3D quasicrystalline materials may exhibit topological 
properties that would have appeared only in 6D periodic 
systems. Furthermore, our approach provides new tools 
for engineering photonic quasicrystals, and especially for 
controlling their surface properties. 



Methods 

The first set of experiments was conducted using an 
array of evanescently coupled waveguides fabricated on a 
semiconductor (ALGaAs) substrate using standard pho- 
tolithography methods. In order to realize the AA model 
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Hamiltonian [cf. equation (1)], we introduced a quasi- 
periodic modulation of the effective refraction index of 
the waveguides while keeping the period constant. The 
average width of the waveguides was 3/im, and the pe- 
riod was 9/im. The sample length was 20mm, which 
is equivalent to 30 tunneling lengths. The light source 
was a continuous-wave laser diode with a wavelength of 
1550nm. The laser beam was focused into a single waveg- 
uide at the input facet of the array using an x40 micro- 
scope objective. The light at the output facet was imaged 
onto an infrared camera (Hamamatsu IR C5840). 



For the adiabatic pumping experiments we used several 
arrays of waveguides written in bulk glas s using femtosec- 
ond laser microfabrication technologySSl, xhe waveguides 
were all identical in both refraction index and width 
2/im), while the inter- waveguide separation was mod- 
ulated in order to realize the "off-diagonal" AA model 
[cf. equation (2)]. The sample length was 75mm, which 
is equivalent to 50 tunneling lengths. The number of 
waveguides in each array was 21. The light source was a 
continuous- wave laser diode with a wavelength of 808nm. 
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SUPPLEMENTARY MATERIAL 



I. THE CHERN NUMBER OF A 
QUASICRYSTAL 

In this section we show that Chern numbers can be 
defined for any given ID quasicrystal (QC). This is ac- 
complished by proving that the Chern density associated 
with a gap of the QC is independent of the phase of the 
quasi-periodicity in the thermodynamic limit. The proof 
is presented for the Aubry- Andre (AA) model [cf. equa- 
tion (1)], but can easily be generalized to other QCs and 
topological indices. 

In order to evaluate the Chern number we consider a 
modified AA model with a similar spectrum. This model 
describes a periodic ID lattice (a ring) of length L with an 
on-site potential parameterized by an irrational number 
b, and the Hamiltonian 

-I- A cos (2TTbLn + (f) tpn , (3) 

where = lb-L\ / L is a rational approximation to b that 
makes the modulation periodic. It is clear that for = 
and L — )• oo this Hamiltonian coincides with the one of 
the AA model. For finite but large L, the spectrum of 
has the same gaps as those of the AA model, up to 
small corrections. 

The phase factor, 0, in the hopping terms introduces a 
phase twist along the ring. Equivalently, it represents a 
magnetic field of 0/27r flux quantum that threads through 
the ring. It plays a role similar to the role of the momen- 
tum in a translation invariant crystal, and is required in 
order for the Chern number to be definecf. 

For an integer 9/21^ this phase factor can be gauged 
out. Hence, the spectrum and the topological properties 
are periodic in 6 modulus 27r. Moreover, since 9 appears 
in the Hamiltonian only via 0/L-terms, perturbation the- 
ory shows that its influence on the spectrum is suppressed 
by a factor of 1 /L. Therefore, for large enough L the gaps 
of -ffft are well defined for every 9. In the thermodynamic 
limit the effect of 9 completely vanishes. 

We can therefore summarize that a gap of TJf, that 
remains open while L is increased corresponds one-to- 
one to a gap of the AA model. Consequently, such a 
gap has, in the thermodynamic limit, the same physical 
properties in both models, in particular the same Chern 
number. 

Consider the shift + 6, where e — 2tiI/L and 
I = 0,l,...,i — 1. There is always some translation of 
the lattice sites n ^ n + rig, where e 0, 1, — 1, 
such that cos (27r6i?i + + e) — cos (27r6i(n 4- n^) + cj)). 
Note that is independent of 0. For a prime L it is 
also guaranteed that there is an corresponding to the 
minimal e = 27r/L. Therefore, if we denote by a trans- 
lation operator by sites, then Hh{cj)+e) = T^Hi,{(j))T~^ . 
The equivalence between the phase shifts and spatial 
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translations implies that the shifts have no physical con- 
sequences, such as closing of the energy gaps. In the 
thermodynamic limit e becomes continuous, which means 
that the spectrum is independent of as mentioned in 
the main text. For finite but large enough L the spec- 
trum weakly depends on (p and the gaps do not close as 
a function of 0. 

A convenient way to evaluate the Chern number asso- 
ciated with a given gap is to consider the projector on 
the states below this gap 



E 



n){n\ 



(4) 



where \n) is an eigenstate oi Hi, with energy i?„, and i?gap 
is the energy at the center of the gap. For large enough 
L, such that the gaps are open for any (p and 0, P{(j), 9) is 
a well defined quantity. Since the projector differs from 
the Hamiltonian only in its eigenvalues, we similarly have 
P(0 -f e) = T^P{(t))T-^. Moreover, by definition 



e) = Um ^ (p(</> + e + A) - P(0 + e)) 



= Hmr,-(F(0 + A)-P(0))T; 



T,d^P{4>)T-^ . 



(5) 



The Chern number associated with the gap is given 
bjl^ 



1 

2TTi 



27r 



d(j)de c{e, (j)) , 



where 



C(0,0)=Tr 



P, 



dP dP 



(6) 



(7) 



is the Chern density. Notably, 



C(0 + e) = Tr ([P(0 + e), [a^P(0 + e), deP{<P + e)]]) 
= Tr (T, [P(0), [d^P{^),deP{4>)]] T,'') 
= . (8) 

We can see that C{(j)) is periodic with 1/L periodicity, 
which means that 



2Tr i-2n/L 



= 27rC(0 = O) + O(l/L). (9) 

This implies that the Chern density, C, is independent 
of (j). Therefore, the integration over in Eq. Q is re- 
dundant, and the Chern number is given by the Chern 
density at a given (j). We can therefore conclude that 
while formally the Chern number is associated with the 
entire family {^^('/>)}0^o' '^^'^ actually associate the 
Chern number with any member of the family H{(f)). 



The generalization of this proof to other QCs is 
straightforward. The proof relies only on the periodicity 
of the Hamiltonian in and on the equivalence between 
phase shifts and translations. Given a Hamiltonian of 
a _D-dimensional QC with d quasi-periodic terms, which 
can be associated with phase shifts, one can define a pe- 
riodic version of the Hamiltonian (similar to Hi,). The 
same conversion from phase shifts to translations still 
holds, as well as the insensitivity to the periodicity and 
phase twists in the thermodynamic limit. Therefore any 
gap of such a I?-dimensional Hamiltonian is character- 
ized by the topological classification of a Hamiltonian of 
D + d dimensions. 

It is worth mentioning that the above argument does 
not imply that the phase (/) has no physical meaning. 
Shifting (/), or equivalently translating the lattice, changes 
considerably the wavefunctions but does not alter the 
energies nor the Chern density. This can be seen by the 
fact that the Chern density is defined as a response to 
infinitesimal changes in 0, but this response does not 
depend on the value of (f). 

As a last note, it can be similarly shown that, in the 
thermodynamic limit the energies and the Chern density 
does not depend on phase twist, 9, as w ell. The proof is 
the same as the one used for the IQHeP^. Therefore, 
the periodic geometry serves only as a formal tool for 
defining the Chern density but, in fact, it is geometry- 
free bulk property. 



II. FURTHER THEORETICAL DISCUSSION 

This section addresses several theoretical issues that 
are mentioned in the main text and deserve additional 
discussion for the benefit of the interested reader. 

A. Appearance of boundary states 

In the main text we state that any given QC, namely 
H{(j)), is associated with a Chern number and belongs 
to a non-trivial topological phase. On the other hand, 
in order to witness the topological gap-traversing modes 
it is necessary to scan (/), i.e. for a given system there 
is either one or no such boundary state on each bound- 
ary. This seems to contradict the conventional notion 
that topological phases, especially those with non-trivial 
Chern numbers, are always accompanied by topologically 
protected gap-traversing boundary states. While this is 
often true, there are examples for topological phases that 
do not necessarily have protected boundary states. 

Consider topological phases that stem from a non-local 
symmetry such as inversion^ or rotation'^. In such sys- 
tems the boundaries may not satisfy the symmetry of 
the bulk, and therefore may not host sub-gap bound- 
ary states. QCs do not rely on symmetry, but on quasi- 
periodicity. Strictly speaking, the boundary breaks the 
quasi-periodicity, hence may not have boundary states. 
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However, the mapping to 2D guarantees that by scanning 
(j) these states must appear. 

Furthermore, in ID topological phases the boundary 
states do not traverse the bulk gap even in the thermo- 
dynamic limit. This is due to the fact that in ID sys- 
tems the boundary does not scale with the system size. 
Therefore at the interface between distinct topological 
phases the number of sub-gap states is of 0(1), e.g. the 
appearance of a Majorana fermion in ID topological su- 
perconductors^. Nevertheless, since our system inherits 
its topological behavior from a 2D system, the boundary 
states indeed traverse the bulk gap, but only when ip is 
scanned. 



B. Quasicrystal vs. crystal 

In our analysis we focused on quasi-periodic systems. 
One may ask what happens to the topological behavior in 
a periodic system. The quasi-periodicity is characterized 
by the irrational number, denoted above by b. This b can 
always be approximated by some rational number p/q, 
where p and q are integers. Such an approximation will 
result in a periodic system. A general property of QCs is 
that this approximation differs from the irrational case 
by 0(1/(7) effects, as we will now demonstrate. 

For the irrational b the spectrum is fractaP^I. Hence, 
it contains gaps of all scales. Conversely, for the rational 
p/q the spectrum is composed of q bandJ^. Therefore, 
taking the rational approximation will cause the gaps 
with energy smaller than (band width)/g to close. 

Moreover, the independence of the energies on 0, that 
is shown in the main text, relies on the equivalence be- 
tween shifts in and translations of the lattice. This 
equivalence holds only for the irrational case. For p/q 
there are only q values of 0-shifts that are equivalent to 
translations. Therefore in the rational case the spectrum 
varies with and is only periodic 27r/g. 

Recalling the proof in Section |Tj it is evident that in 
the rational case also the Chern density depends on (p 
with 27r/g periodicity. Therefore, in order to evaluate 
the Chern number one has to perform the integration 
over (j). This means that the Chern number is associated 
only with the whole family of ID systems, and that a 
single periodic system belongs to the trivial phase. This 
can be also understood by the fact that deformation of 
such a family into a trivial family will cause the bulk gap 
to close only for q values of 0, while for the other values it 
remain open. Nevertheless, in case that the whole fam- 
ily of Hamiltonians is non-trivial, boundary states will 
indeed traverse the bulk gap as a function of 4>. 

Finally, we note that for a finite system of length L, 
it will be physically impossible to distinguish between an 
irrational b and a rational approximation p/q for q^ L. 



C. Robustness 

Topological phases are characterized by having their 
properties being robust to perturbations. This is also 
the case for QCs. The topological classification of QCs 
holds also in the presence of disorder, as long as the disor- 
der preserves the quasi-periodicity and does not close the 
bulk gap. Therefore (i) continuous deformation between 
two topologically distinct systems results in a phase tran- 
sition also in disordered systems; (ii) The appearance of 
boundary states that traverse the gap as a function of 4> 
is also guaranteed. Since, any experimental setup is in- 
herently disordered to some level, this property was also 
demonstrated in our experiments. 

Considering ferniionic systems, temperature should 
also be taken into account. Since temperature only deter- 
mines the occupation, the topological phase is stable to 
temperatures that are smaller than the bulk gap. Simi- 
larly, when interactions are introduced a quantized Chern 
number is still expected to hold as long as the bulk gap 
remains open. 



D. Alternative gedanken experiment 

The association of a Chern number with a given QC, 
together with the above discussions, can be demonstrated 
in the following gedanken experiment. Suppose we have 
a sample of a ID QC of finite but large length that is pop- 
ulated with electrons up to some chemical potential. The 
sample should be put under an STM, which measures its 
local density of states (LDOS). In general, we will find 
that the spectrum is gapped over the bulk, whereas at 
each boundary a sub-gap boundary state may or may not 
appear. If there are no boundary states, can we conclude 
that the QC is trivial? 

The answer is no. We should cut the very last piece 
of the lattice at one of the boundaries, and repeat the 
measurement. The gap at the bulk will, of course, remain 
the same, but a sub-gap boundary state will probably 
appear. We repeat this process several times. If the 
state is topological, then during repeated cuts the state 
will disappear and reappear, but each time at a different 
energy. This happens since cutting of the last piece can 
be thought of as a shift in (j). 

Integrating the LDOS at the boundary measured dur- 
ing this process, the gap will be gradually filled. Notably, 
given the quasi-periodicity 6, each cut is equivalent to a 
shift of 27Tb. Therefore, by restoring the sub-gap spec- 
trum as a function of cf) the Chern number can also be 
resolved from the measured data. 

At this point, one should notice that in case of a pe- 
riodic system with frequency p/q, boundary states will 
also disappear and reappear, but now their energy will 
have only q different values. Hence, the integrated LDOS 
within the gap will be composed of q peaks, rather than 
being continuous. 
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The described experiment promotes several observa- 
tions: First, It is performed on a single member of the 
QC family with open boundaries. We can therefore see 
that it is a bulk property and that the periodic geometry 
that was introduced in Section [l] is merely a formal tool. 
Second, if we take another member of the family (i.e. a 
sample with the same b but a different 0) at the same 
chemical potential and use this protocol, we would ex- 
tract the same Chern number. This demonstrates the 
independence of the Chern density on (j>. Third, the 
boundary states that fill the gap in the integrated LDOS 
originate from topological edge states in 2D. Last, it is 
evident that the physics is robust to disorder and tem- 
perature that does not close the bulk energy gap. 



such QCs. This is done by defining a smooth version of 
the projection procedure. This procedure keeps the en- 
ergy gaps unchanged, and therefore has the same topo- 
logical properties. For example, in ID, a QC can be ob- 
tained by projecting a square lattice on the line y = x/t, 
where r is irrational. The normalized spacings between 
the lattice sites n and n -|- 1 is therefore 



1 + (r - 1) 



(n + 2) 



(n + l) 



(10) 



The case of r = (l + v5) /2 is the well studied Fibonacci 
qCSIIsi smoothed projection can be of the form 



III. GENERALIZATION TO OTHER TYPES OF 
QUASICRYSTALS 

In the main body of the paper we considered two exam- 
ples of ID quasicrystals that exhibit topological proper- 
ties: (i) the "diagonal" Aubry-Andre model [cf. equation 

(1) ] in which the quasi-periodicity is in the on-site po- 
tential, and (ii) its "off-diagonal" version [cf. equation 

(2) ] that has its quasi-periodicity in the hopping term. 
We then argued that tight-binding Hamiltonians of other 
QCs and of higher dimensions also have such topological 
properties. 

However, there is a cla ss of QCs that is obtained by a 
"cut and project" methocPE^, for which one might think 
that our argument does not hold. For these QCs the 
additional d.o.f. are shifts of the "cut" in the higher- 
dimensional space. Such shifts modify the Hamiltonian 
in a discontinuous manner. Hence, as our argument relied 
on the continuity of the on-site and hopping terms as a 
function of the additional d.o.f. , it seems that these QCs 
do not fit our framework. 

We now show that our generalization indeed includes 



1 — tanh/3 



1 



(11) 





'2n + 3 1' 








+ 01 + COS - 


) 




2r 4 







. dn (/3, = 0)= d„. 
in a way that is con- 



It can be easily shown that limj_ 
Note that is embedded into dn, 
sistent with our general argument. 

The spacing between the lattice sites enters the Hamil- 
tonian through the hopping terms. Since a Hamiltonian 
with dn can be continuously deformed into one with dn 
without closing the energy gaps, they belong to the same 
topological class. Therefore, given a QC which is ob- 
tained by a "cut and project" procedure, one should find 
a smoothed version of the projection and accordingly ob- 
tain a smoothed Hamiltonian. Now, the topological prop- 
erties of the "cut" Hamiltonian can be deduced from the 
smoothed one. Note, that for "cut and project" quasi- 
periodicity which is incorporated in the on-site terms, 
a similar smoothing can be performed. Thus, we have 
shown that the "cut and project" QCs also fit into our 
general argument. 
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